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We present the LHC predictions for the WHDG model of radiative, elastic, and path 
length fluctuating energy loss. We find the p± dependence of Raa is qualitatively very 
different from AWS-based energy loss extrapolations to the LHC; the large reach of 
the year one data at the LHC should suffice to distinguish between the two. We also 
discuss the importance of requiring a first elastic scatter before any medium-induced 
elastic or radiative loss occurs, a necessary physical effect not considered in any previous 
models. 

1. Introduction 

Currently, several claim that their respective descriptions of en- 

ergy loss properly account for the high-p^ suppression of central pion data observed 
Y Q q 

at RHIC ° Other RHIC observables have been suggested as a means of differ- 
entiating between these, such as their predictions for the centrality dependence, 
back-to-back jet quenching, two and three particle correlations, photon spectra, 
etc. We suggest that the large p± lever arm of pion data at the LHC wil l distin- 
guish between the WHDG and the AWS models. We will not discuss further; 
higher twist predictions for the LHC do not seem to have been calculated, and 
AMY LHC results go out to only low p±_ 30 GeV). 

The important quantity in any of these calculations is P(e), the probability that 
the final momentum is some fraction of the initial momentum, pj_j = (1 — £)p±,i- 
For calculations that include only radiative processes, P{e) — Pradi^), where 

j P^d{x)+Prad{x) 
I n 

The possibility that no gluons are emitted (subsequent to the original gluon radi- 
ation created by the initial hard scatter) is encapsulated in the coefficient, Pq , of 



9, 2008 0:55 WSPC/INSTRUCTION FILE horowitz 



2 William A, Horowitz 

the S{x) term. 

Mustafa showed that the elastic contribution to charm energy loss is comparable 
to the radiative losJl^. This work was extendeci^ to show that, despite long held 
assumptions to the contrary, elastic energy loss is in fact not negligible for all parton 
jet flavors. The WHDG model includes elastic energy losses in its Raa calculations 
by convolving the two probability distributions, 

P(e) = j dxPrad{x)Pei{e-x). (2) 

However, none of these previous models include the possibility that the parton 
jet escapes the medium completely unmodified, Pq = exp(— A^c), where Nc is the 
average number of elastic collisions suffered along the parton path. Note that this 
is different from, and in addition to, the probability of radiating no gluons, Pq — 
exp{—Ng), included in the Prad term in Eq. ([T]). 

Due to the approximations used in the radiative calculations overabsorption, 
P(e > 1), has a large support for highly suppressed jets. Usually, one of two pre- 
scriptions is applied to remove this unphysical artifact. Either the integrated excess 
probability weighs an explicit delta function, 

/oc 
dxPoid{x)6{l~e), (3) 

or reweighs (rw) the original distribution, 

dxPoid{x) 

Clearly the latter approach leads to larger Raa values for the two. For large over- 
absorption energy loss details are lost to the removal of the unphysical e > 1 region. 



2. WHDG Model 

In our extended theory of elastic loss in addition to radiative loss calculated in a 
realistic geometry, we take 

Raa{p±) « ( / de{l - e)"(f-)p(e)),eo™. (5) 



n{p±) comes from the power law approximation to the pQCD production spec- 
trum (minus either one or two depending on the Jacobian). P(e) is the convolved 
probability distribution P(e|x, (/i>,p^) = / dxPradix\x,(l),p±)Pci{e — x\x,(j),p±). And 
(■ ■ ■)geom corrcsponds to geometrical averaging. 



_ JdL...SiL-L) 



see^ and references therein for details. 
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The first line of Eq. ([S]) is a simplification of full geometrical averaging, and 
presents one with a number of choices for L. The three most natural possibilities 
are 

{Lprod 
Lobs{p±) 
Lftt, 

where Lprod is the average length that all produced hard partons travel, Lprod = 
^ JcPxd(j)pprodL{x,(j)), Lobs{p±) is the average length that an observed jet trav- 
elled on its way out of the medium, Lobs{p±) = ^ J(fixd(j)pprodL{x,(j)) /de(l — 
e)^(p±)P{e\x, 4>,p±), and Lfa is the single fixed length that best approximates the 
full geometry average, given by the second line of Eq. While it is generally im- 
practical to use Lfit in calculations (one must first compute the proper geometry 
average first anyway in order to find L/it), one sees from Fig. [l]that employing a 
single length without calculating the full geometry average is problematic; there is 
no a priori method to estimate the degree to which using a single, representative 
length deviates from the full solution. In fact, Fig.[l]shows the Lfu hierarchy differs 
from one's naive expectations. 



(7) 




Fig. 1. All plots depict midrapidity for most-central collisions, (a) The histogram gives the dis- 
tribution of path lengths traversed by hard scatterers. The lengths, L(a;^,(/>), are weighted by 
the probability of production and averaged over azimuth; these quantities were found using the 
Glauber thickness function from a realistic Woods-Saxon base geometry. A purely geometric quan- 
tity, P{L) is the same for all jet varieties. However, no single Lfit can best reproduce the results of 
the full distribution for all parton jet flavors, (b) The WHDG model reproduces the PHENIX pion 
Raa{p±) data^'®'^ for the conservative, entropy-constrained input parameters of dNg/dy = 1000 
and as = .3. Model sensitivity to varying the medium density is much greater than that observed 
in ^; however, the large sensitivity to changes in ots suggests a need for relaxing the fixed coupling 
constant approximation in future work. 
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3. RHIC Results 

We see from Fig. [T] that for the conservative as of .3 and entropy-constrained 
dNg/dy = 1 000, t he WHDG model is consistent with the RHIC pion data observed 
by PHENIXPI. 

While saw these results as sensitive enough to changes in the 
input to be considered not fragile (unlike the conclusions of ^ , a more sophisticated 
statistical analysis showed that, when requiring the rather high 90% confidence 
level, the combination of theoretical sensitivity and experimental error results in 
an approximately factor of 2 range of allowable medium densitiei^^l, We note that 
previous, radiative only GLV-based energy loss models also are consistent with the 
datal^ in comparison, the addition of elastic losses does not oversuppress our results 
for two reasons: we also include geometrical path length fluctuations, and we use 
a fixed as - ^ incorporates some running of the coupling. Due to the low q± and 
k±'s involved, a running will often be evaluated at its cutoff; for SI, a^"'^ — .5. 
We note that the inclusion of elastic loss allows for a small Raa without much 
overabsorption in Prad, Pel, or their convolution. 

In order to understand the qualitative differences that will arise i n the LHC pre- 
dictions, we briefly discuss the RHIC results from AWS-type modelJ2E!_ Xo fit the 
pion suppression, a nonperturbatively large q is used. Specifically, q is proportional 
to the three-fou rth 's root of the energy density of the med ium , q = cc"!^ ^ where 
pQCD estimateJi^ give 2\ rather, c ~ 8 — 20 is needed in 1^121 This is due to the 
combination of radiative only energy loss, oversirnplified treatment of geometrjUJJ, 
and the unitarity-violating Pq > 1 in their model^. 



4. LHC Predictions 

We show in Fig. [2] the large qualitative difference in dependence for predicted 
LHC pion Raa from the WHDG model and from two different implementations 
of the AWS model; one easily sees the dramatic rise in Kaa with increasing 
from the first as opposed to the flat in pj_ results of the latter two. The consistency 
of Vitev's curve^ wit h ou rs over a range of dNg/dy and the consistency of the 
two AWS calculation^2l3| suggest that this is a robust result. The origin of the 
difference can be easily understood given our discussion of RHIC results. For the 
case of WHDG, the pion Raa at RHIC does not require much overabsorption. The 
modest increase in me dium density, 2 — 3 based on eith er an extrapolation of 
PHOBOS resultJ14il5| 

or predictions from the for the LHC leads to 

small energy losses at high momenta that can be well approximated by the pocket 
asymptotic energy loss formulae 

erad = AErad/E ^ \og{E / L) / E (8) 
tei = AEei/E - \og{VET/mg)/E. (9) 

As p± increases the \og{E)/E reduction in energy loss is not compensated by the 
slow (in comparison to RHIC) increase in the power law, n{p±), partonic production 
spectrum; thus Raa increases with p±. On the other hand, the AWS mimics the 
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small normalization of the RHIC data by highly suppressing their jets; the details 
of energy loss are lost (in the delta functions at zero and one, or at zero and the 
reweighting of P), thus flattening the results. Moreover, the two AWS models rep- 
resented in Fig. [2] used EKRT-type medium density scalinjl^; this makes the LHC 
~ 7 times more dense than RHIC. The LHC jets are thus even more dramatically 
overabsorbed and, again, flat in p± . 
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Fig. 2. LHC predictions for several energy loss models. WHDG curves correspond to ^, GLV 
to PQM AWS to 3, and AWS to 2. For the latter two, rw indicates the use of Eq. ^ for 
overabsorption; otherwise Eq. (|3]l was employed. Also for the latter two, q = 100 and q = 68 were 
the values of the AWS input parameter, respectively. Notice the sharp rise in p± for the WHDG 
and GLV curves as opposed to the flatness of the AWS and PQM AWS results. 

5. Probability of Nothing 

An outgoing parton jet must encounter at least one scatter in order for medium- 
induced energy loss to occur. There exists, then, a new (in the sense that it has not 
been considered before) probability of no in-medium energy loss, Pq — exp(— iVc), 
where N^. = ^ dza ei{z) p{z) is the average number of elastic collisions for the jet. 
Note that Nc has a strong dependence on the Casimir associated with the quark or 
gluon jet. After this first elastic collision the total probability of energy loss is then 
the convolution of the elastic and radiative processes: 

F(e) = Pa5{e) + (1 - Pq) j dxPrad{x)Pei{e - x). (10) 

In the above equation, this Po5{e) piece is in addition to the probability of no glue 
emission, P^5{x), in Prad{x), Eq. ([1]). 

For fixed as = .3, including Pq physics accounts for 50% of Raa- Allowing as{T) 
to run as asiq^ = 2ttT{z)) reduces Pq by a factor of 2. Finally, integration over 
momentum transfers with as{q'^) given by vacuum running formally gives Pq — 0. 



WHDG dNj/dy = 1750 

WHDG dN /dy = 2900 



- GLV dNj/dy = 2000 

GLV dNj^/dy = 3000 




± 
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6. Conclusions 

We see from Fig.[2]that the LHC Raa{p±) pion data will distinguish between energy 
loss models. While WHDG, with its radiative, elastic, and path length fluctuations, 
predicts a significant rise as a fmiction of AWS-type models predict flat p±_ 
dependence. This flatness in pj_ of the latter is due to the high suppression of their 
jets, which forces the energy loss details into the region of overabsorption. Contrar- 
ily, the former predicts a much smaller fractional energy loss for jets; the high-px 
behavior is well approximated by the analytic asymptotic energy loss formulae and 
is thus responsive to their details. We also found that moderate in opacity Raa 
predictions are sensitive to noninteracting free jets, whose influence has not been 
considered previously. 
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